Double-exchange is not the cause of ferromagnetism in doped manganites 
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The coexistence of ferromagnetism and metalhc conduction in doped manganites has long been 
explained by a double-exchange model in which the ferromagnetic exchange arises from the carrier 
hopping. We evaluate the zero-temperature spin stiffness -D(O) and the Curie temperature Tc on 
the basis of the double-exchange model using the measured values of the bare bandwidth W and 
the Hund's rule coupling Jh- The calculated D(Q) and Tc values are too small compared with the 
observed ones even in the absence of interactions. A realistic onsite interorbital Coulomb repulsion 
can reduce -D(O) substantially in the case of a 2-orbital model. Furthermore, experiment shows that 
-D(O) is simply proportional to x in Lai-ajSrjjMnOa system, independent of whether the ground 
state is a ferromagnetic insulator or metal. These results strongly suggest that the ferromagnetism 
in manganites does not originate from the double-exchange interaction. On the other hand, an 
alternative model based on the d — p exchange can semi-quantitatively explain the ferromagnetism 
of doped manganites at low temperatures. 



The discovery of "colossal" magnetoresistance in thin 
films of the manganite perovskites Rei_a;Da;Mn03 (Re — 
a rare-earth element, and D = a divalent element) [QJ|] 
has stimulated extensive studies of magnetic, structural 
and transport properties of these materials js). The co- 
existence of ferromagnetism and metallic conduction has 
long been explained by the double-exchange (DE) model 
where the effective hopping for the manganese 3d 
conduction electrons varies with the angle between the 
manganese core electrons due to a strong Hund's cou- 
pling. However, Millis et al. ||] proposed that, in addition 
to the double-exchange, a strong electron-phonon inter- 
action arising from a strong Jahn- Teller effect should be 
involved to explain the basic physics of manganites. In 
this modified model, the primary cause of the ferromag- 
netism of doped manganites is still the double-exchange 
interaction. 

In the DE model, it is implicitly assumed that doped 
carriers are Mn Cg electrons. This assumption is not jus- 
tified by both electron-energy-loss and photoemission 
spectroscopies which have shown that the ferromag- 
netic manganites (a;<0.4) are doped charge-transfer in- 
sulators with carriers mainly residing on the oxygen or- 
bitals. Now a question arises: Does the ferromagnetism 
of doped manganites really originate from the DE inter- 
action? If not, what causes the ferromagnetism in these 
compounds? One way to address this fundamental issue 
is to make a quantitative comparison between the pre- 
dicted properties of the DE model and experiment. 

There are two important parameters in the DE model, 
namely, the bare bandwidth W of the Cg bands, and the 
Hund's rule coupling Jh between Cg and t2g electrons. 
These parameters are related to an optical transition be- 
tween the exchange splitted Cg bands and thus 
can be determined from optical data. With these unbi- 



ased parameters, one can calculate the zero-temperature 
spin stiffness -D(O) and Tq- Here D{Q) is defined as 
ujq ~ D{0)q^ with ujq being the magnon frequency. When 
one introduces interactions such as electron-phonon and 
electron-electron interactions, the magnitudes of both 
D{0) and Tc are generally reduced. 

Here we use the measured values of the bare band- 
width W and the Hund's rule coupling Jh to calcu- 
late £'(0) and Tq on the basis of the double-exchange 
model. The calculated 1^(0) and Tc values are too small 
compared with the observed ones. A realistic onsite in- 
terorbital Coulomb repulsion can reduce D{Q) substan- 
tially in the case of a 2-orbital model. Moreover, exper- 
imental data show that D{0) is simply proportional to 
X in Lai_a;Srj;Mn03 system, independent of whether the 
ground state is a ferromagnetic insulator or metal. These 
results provide strong evidence that the ferromagnetism 
in manganites is not caused by the DE interaction. On 
the other hand, an alternative model based on the d — p 
exchange can well explain the ferromagnetism of doped 
manganites. 

Now we start with a Kondo-lattice type Hamiltonian 
pl| , which leads to Zener's DE model when Jh — > oo, 

<ij>aba 

—Jh ^ Si ■ dl^^cjapdiap + Hint- (1) 

iaaf3 

Here df^^ creates an electron in eg orbital a with spin a, 
tlij is the direction-dependent amplitude for an electron 
to hop from orbital a to orbital 6 on a neigboring site, 
and Hint represents the other interactions. The calcu- 
lated band structure is well fit by a tfj*, which involves 
only nearest-neighbor hopping that is only nonzero for 
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one particular linear combination of orbitals, i.e., t^^ cx 
t, where t is a characteristic hopping ampUtude that is 
related to the bare bandwidth W hy W = At. Here we 
still call Eq. 1 as 2-orbital DE model rather than 2-orbital 
Kondo-lattice model for convienence. The quantum and 
thermal average of the hopping term in Eq. 1, defines a 
quantity K: 



K = (l/6iV,,te) 



<ij>aba 



ttHdLd,bo. + h.c.), (2) 



The quantity K is related to the optical spectral weight 
by a familar sum rule, 



K = [ duj(Ti{uj), 

Jo 



(3) 



where CTi is the real-part optical conductivity contributed 
only from the Cg electrons, and Oo is the lattice constant. 
The quantity K generally consists of the Drude part Kjj, 
and incoherent part Ki which, in general, involves inter- 
band and intraband optical transitions. The Drude part 
Kd can be related to the plasma frequency Qp as 



Kd 
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On the basis of Eq. 1, Quijada et al. [|Tl[ showed that, 
to the order of \/ Jhi the spin stiffness D{Q) is, 



45* ^ 



JhSK' 



(5) 



where S = 3/2, S* ^ S + {1 - x)/2, and r; = 1.04 when 
Hint = 0. The presence of interactions may change the 
value of rj. A similar result was obtained by Furukawa 
1^ for an 1-orbital DE model using the dynamical mean 
field method, but the value of 77 is doping dependent and 
less than 1. 

From the above equations, one can calculate K and 
D{Q) using realistic values of the bare bandwidth W and 
the Hund's rule coupling Jh- Both the local density 
approximation (LDA) and "constrained" LDA 
calculations show that — 1.5 eV. The calculated Jh 
value is very close to the atomic values for 3d atoms. This 
is reasonable because Jh is not screened when the ion is 
put in a solid. The bare bandwidth W cannot be cal- 
culated reliably for 3d-metal based compounds due to a 
strong correlation effect. Fortunately, the values of both 
W and Jh can be determined from an optical transition 
between the exchange splitted eg bands |^,|o| . The peak 
position of this optical transition is about 2Jh, and the 
width of the peak contains information about the bare 
bandwidth ||,|lO|. From the optical data of Ref. ||ll|,|l3l, 
one finds Jh — 1.6 eV and W = 1.6-1.8 eV by compar- 
ing the data with the theoretical predictions 1^,^ . The 
value of Jh obtained from the optical data is in excellent 
agreement with the calculated one. This implies that the 



feature appeared at about 3 eV in the optical data indeed 
arises from the optical transition between the exchange 
splitted Cg bands. 

The quantity K° for noninteracting 2-orbital model 
can be evaluated when the bare bandwidth W is known. 
Takahashi and Shiba have calculated the optical con- 
ductivity using a tight binding (TB) approximation of 
the band structure. From their calculated result for the 
interband optical conductivity, we evaluate that Kj ~ 
0.088t for X ~ 0.3. One should also note that the mag- 
nitude of t defined in Ref. is 1.5 times smaller than 
the t defined here. Since = l.2K° [|l5|, then K^, = 
0.106t and K° = 0.194^. The LDA calculation for a cubic 
and undistorted structure shows that |l^ W = 3 eV and 
hn° = 1.9 eV. Using Eq. 4 and hn° = 1.9 eV, one yields 
= 78.6 meV. Since i = Vl^/4 = 0.75 eV, one readily 
finds that K"^ = 0.105t, in remarkably good agreement 
with that {K"^ = 0.106i) estimated from the TB approx- 
imation. This justifies the relation K° = 0.194i obtained 
from the TB approximation. It is interesting to compare 
the present result with those reported in Ref. ||l^ and 
111. In Ref. 0, it is found that K° = 0.34i for x = 0.3 
using the dynamic mean field method. In Ref. ||ll[| , Qui- 
jada et al., claimed K° — 0.46i for x — 0.3, which might 
be true if the two eg bands have the same dispersion. 
Therefore the quantity K° is significantly overestimated 
in Ref. 0] and |ll|. 

When the Hund's coupling Jh is turned on, the quan- 
tity K is reduced compared with K° . For Jh = 00 and 
X = 0.3, K = 0.77 K° |lO|jl^. It was also shown that 0] 
the reduction factor (0.77) is basically the same for Jh 
> t. Therefore, we have K = 0.147i for x = 0.3. Us- 
ing t = 0.4 eV, Jh = 1.6 eV, and K = 0.147i, we yield 
D{0) = — 25 mV from Eq. 5. The negative value of 
£'(0) implies that the ferromagnetism is not sustainable 
with these unbiased parameters. It might be possible to 
have a small positive -0(0) if one includes higher order 
terms in Eq. 5. Nevertheless, the theoretical D{0) value 
is too small compared with the measured ones (160-190 
mV A2) |l|-|§. 

Now we turn to the calculation of Tc for x — 0.3. For 
an 1-orbital DE model with Jh = 00, the dynamic mean 
field (DMF) calculation shows that T*^^ = 0.078t/A:Bji, 
while Monte Carlo simulations yield Tc — O.Oit/kB pi |, 
or Tc = O.OM/kB This implies that the DMF calcu- 
lation overestimates Tc by a factor of about 2 due to the 
neglection of fiuctuations. We would like to mention that 
the magnitude of t defined in Ref. |2^ and is 3 times 
smaller than the t defined here. For realistic parameters t 
= 0.4 eV, = 4i = 1.6 eV, the DMF calculation obtains 
T^'F = 0.038t/A:B = 180 K |]. Considering the fact that 
the DMF method can overestimate Tc by a factor of 2, 
one has Tc ~ 100 K. For the 2-orbital DE model, the 
DMF calculation shows T*^^ = 0.07i/fcs = 324 K. 
Since the quantity K calculated in Ref. pOl is overesti- 
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mated by a factor of about 2 (as discussed above), and 
the DMF method itself can overestimate Tc by about 2 
times, the real Tc should be about 100 K, which is com- 
parable with the value for the 1-orbital model. This is 
reasonable because the bare values of K° for both 1- and 
2-orbital models happen to be very similar in the case of 
X = 0.3. Therefore, both the 1- and 2-orbital DE models 
cannot explain the observed Tc values with the unbaised 
parameters. 

The above calculations have not taken into account any 
other interactions such as electron-phonon interaction 
and electron-electron correlation. The electron-phonon 
interaction can substantially reduce the K and thus D{Q) 
if the coupling constant A = Ep/2t is greatly larger than 
1 | p3| , p^ . In reality, the polaron binding energy Ep in 
manganites is estimated to be about 1 eV So A ~ 
1, which suggests that the electron-phonon couphng does 
not lead to a sizable decrease in 15(0). On the other hand, 
the electron-electron correlation in the 2-orbital model 
can lead to a large reduction in K , as demonstrated by 
Horsch and coworkers . The parameter U = U' - Jab 
(where U' is the onsite interorbital Coulomb repulsion, 
and Jab is the interorbital Hund's coupling) has a strong 
influence on the value of K . With realistic values of 
?7 = 3 eV |6|,|3 and t = 0.4 eV, K can be reduced by 
a factor of about 2 compared with the bare one . Al- 
though their results are for 2-dimensional finite clusters, 
we would expect a similar result for a real 3-dimensional 
system. 
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FIG. 1. Zero-temperature spin stiffness -D(O) vs x in 
Lai-a^Sr^^MnOa system. The -D(O) value for x — 0.15 is taken 
from Ref. [^, and the other data from Ref. |2^]. It is evident 
that -D(O) is simply proportional to x, independent of whether 
the ground state is a ferromagnetic insulator (e.g., x = 0.15), 
or a ferromagnetic metal (e.g., x = 0.175, 0.20, 0.30). 

In addition, it is striking that the observed spin stiff- 
ness varies only with x, namely, D(0) cx x, as seen clearly 
from Fig. 1. The doping dependence of -D(O) shown in 
Fig. 1 could be qualitatively explained by the 1-orbital 
DE model §,|§ where D{0) (x K x{l-x) for Jh ^ oo 



pO| , as plotted in Fig. 2. Nevertheless, the doping depen- 
dence of D(0) is very different from that prediced from 
the 2-orbital model. By analogy to the 1-orbital model, 
one can easily show that D{0) (x K oc {1 — x){l + x) for 
the 2-orbital model, as demonstrated in Fig. 2. Here we 
have implicitly assumed that the two Cg bands have the 
same dispersion. The solid circles in Fig. 2 are the results 
evaluated using the dynamical mean field method with 
Jh — oo |l^. It is evident that the analytical expres- 
sion is in excellent agreement with the numerical result. 
Comparing the results shown in Fig. 1 and Fig. 2, one 
clearly sees that the 2-orbital model cannot explain the 
observed doping dependence of -D(O). 




FIG. 2. The doping dependence of -D(O) predicted by the 
1-orbital and 2-orbital DE models with Jh —> oo. The solid 
lines represent -D(O) oc x{l — x) for the 1-orbital model, and 
D{0) oc (1 — a::)(l + 3;) for the 2-orbital model. The solid circles 
are the values numerically calculated from the 2-orbital model 




The question is why the DE model cannot explain 
the ferromagnetism in doped manganites. As mentioned 
above, the DE model implicitly assumes that doped car- 
riers are Mn Cg electrons, which is not the case accord- 
ing to the electron-energy-loss and photoemission spec- 
tra |7|Jq| . Furthermore, the "constrained" LDA calcula- 
tion [n3[ shows a large onsite Coulomb repulsion of about 
8-10 eV, in agreement with the photoemission data 
The simple LDA calculation which ignores the strong cor- 
relation effect shows a bare plasma frequency of about 
1.3 eV for x — 0.33 with a distorted structure deter- 
mined by neutron scattering pi] . The bare plasma fre- 
quency calculated is much smaller than the one observed 
in Ndo.TSro.aMnOa (3.3 eV) [|2). The large bare plasma 
frequency observed in this material is consistent with the 
fact that doped holes reside mainly on the oxygen orbitals 
with a large bandwidth. The bare plasma frequency of 
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about 3.3 eV for single conduction band (oxygen band) 
implies a bare K° of about 0.24 eV, which gives an up- 
per limit for the K in the presence of interactions. The 
electron-phonon interaction with a coupling constant A 
~ 1 will reduce the K slightly 24|, but can signifi- 
cantly decrease the Drude weight. Optical data indeed 
show that the K for Ndo.rSro.sMnOs is about 0.2 eV 
[ pl| , while the effective plasma frequency is about 0.57 
eV Q. The Drude weight is reduced by a factor of 33, 
implying small polaronic carriers in the low-temperature 
ferromagnetic state. 



A = 4.5 eV 




FIG. 3. The schematic band structure for doped man- 
ganites (a::<0.4) constructed from the LDA -I- U calculation 
[ L3 . The energy scales are consistent with the optical data 
[ 11 ]3^ , ^ . The density of the oxygen holes is equal to x per 
cell while the densities of the Mn'^"'' and Mn*^ ions are the 
same due to the charge disproportion (2Mn'^^ — > Mn^+ + 
Mn''+) le). 

In Fig. 3, we plot a schematic band structure for doped 
manganites {x< 0.4), which is extracted from the LDA 
-|- U calculation. Here we have assumed that the local 
Jahn- Teller distortions still survive upon doping, but the 
average magnitude of the distortions decreases, in agree- 
ment with the experiments The doping with a 
divalent element shifts down the band due to the de- 
crease of the Jahn- Teller distortions. The density of the 
oxygen holes is equal to x per cell, while the electron 
carrier density in the majority band (corresponding 
to the density of the Mn^+ ions) is the same as the hole 
carrier density in the majority band (corresponding to 
the density of the Mn''^+ ions). In other words, the dop- 
ing does not change the average valence of the Mn ions, 
but lead to the charge disproportion (2Mn'^+ Mn^+ -|- 
Mn'*+). This is because the quench of the static Jahn- 
Teller distortions by doping makes the Mn^+-Mn^+ pairs 
more stable than the Mn3+-Mn3+ pairs [ ^6| . The current 
band structure is consistent with the optical transitions 
at the photon energies of about 1.5 eV, 3.0 eV and 4.5 
eV (ll]||0. The optical spectral weight for the 4.5 



eV transition should be much larger than for the 1.5 eV 
transition, as observed js^ . This is because the unoccu- 
pied state density for the former optical transition (i.e., 
3 minority t2g and 2 minority tg states per cell) is at 
least 5 times larger than that for the latter one (i.e., less 
than 1 majority eg states per cell). The optical transi- 
tion at about 3 eV is related to the transition between 
the exchange splitted Cg bands. 

What is an alternative model for the ferromagnetism in 
doped manganites? If we consider an oxygen hole (spin 
1/2) sitting in between two Mn ions, an exchange inter- 
action between the oxygen and Mn spins {d—p exchange) 
will lead to a ferromagetic interaction between Mn spins 
p8| . In this case, the ferromagnetic exchange energy be- 
tween two Mn spins is J = C(ip^/A^(), where C is a 
numerical factor of order of 1, tpd is a hybridization ma- 
trix element between the d and p orbitals, and Act is a 
charge transfer gap. The long-wave spin stiffness for a 
cubic material is (H] L»(0) = {S* Ji3\Ri - 
where Jy is the exchange energy between pairs of spins 
at sites Ri and Rj, and the sum is over six neighbor- 
ing sites. Clearly, J if there is an oxygen hole in 
between two Mn ions, while = if not. Since the 
oxygen hole density is equal to x per cell, then D{0) = 
{2/3)xJS*al, independent of whether the holes are lo- 
calized or not. This mechanism can naturally explain 
why £'(0) is simply proportional to x, as shown in Fig. 1. 
Using the Slater-Koster parameter Vpda = 1-8 eV ob- 
tained from the analysis of photoemission data [||, or 
from the TB fit to the LDA band structure , we yield 
tpd - -iV^/2)Vpda -= -1.56 eV g. Taking A^t = 4.5 
eV [|,0|37|, tpd =-1.56 eV, C = 1, and a; = 0.3, we ob- 
tain D{0) — 0.5 eV A^, which is larger than the observed 
values by about three times. This is not unreasonable 
since the numerical factor C might be less than 1. A 
more quantitative calculation of J based on the d ~ p 
exchange is essential to address this issue. 

It should be noted that the ferromagnetic "bonds" with 
the exchange energy J will be randomly distributed over 
the real space if the oxygen holes are localized, whereas 
the distribution of the bonds will become more homoge- 
neous when the holes are more mobile (a motion narrow- 
ing effect). This can account for a conventional spin- wave 
dispersion in high Tc materials ]39| ] where the conductiv- 
ity is high, and an unconventional magnon softening near 
the zone boundary for lower Tc compounds [ p^ where 
the conductivity is lower. When the oxygen holes are or- 
dered (charge-ordering), as observed in Lao.ssSro.isMnOa 
p8[ , the ferromagnetic bonds will form a "superlattice" , 
leading to a splitting of spin-wave dispersion at a wave 
vector equal to the incommensurability of the superlat- 
tice as observed p8| . 

The above simple model would imply that the Curie 
temperature Tc should be proportional to x. In reality, 
Tc is strongly dependent on the cation radius va of the 
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Tiei-xT)x site of the perovskite Rei_2;Da;Mn03 even if x is 
fixed Moreover, a giant oxygen isotope shift of Tc 
has been observed in these compounds 44|. There- 
fore, the simple d — p exchange model cannot account 
for these unusual phenomena. Recently, Alexandrov and 
Bratkovsky have proposed that, in addition to the 
d — p exchange interaction, there is a strong electron- 
phonon interaction that may lead to the formation of 
small (bi)polarons. In the paramagnetic state, the singlet 
bipolarons (spin zero) are stable and the ferromagnetic 
interaction is produced by the thermally excited polarons 
(spin 1/2). Thus, the Curie temperature Tc is selfconsis- 
tently determined by the polaron density at Tc- Within 
this scenario, Tc can strongly depend on the electron- 
phonon coupling strength and the isotope mass | |45[ |, in 
agreement with experiment 44 1. At zero temperature, 
all the carriers will be polarons if JpdS is much larger 
than the bipolaron binding energy A, where Jpd is the 
exchange energy between Mn and oxygen-hole spins. In 
this case, D[0) oc x, as in Lai-^^Sr^-MnOa system. On the 
other hand, there is a mixture of polarons and bipolarons 
if JpdS is slightly larger than A. Then the relation D{0) 
oc X does not hold, and D{0) is proportional to the den- 
sity of polarons which is less than x per cell. One should 
also expect that the mixture of mobile polarons and lo- 
calized bipolarons would lead to a dynamic/static phase 
separation since bosonic and fermionic carriers tend to 
separate in real space. 

Finally we would like to address whether the band 
structure shown in Fig. 3 can be consistent with an effec- 
tive 1-orbital DE model with a large t ~ 1.2 eV. We can 
show that in order to have a Tc of 380 K with such a t, 
one needs an unphysically large value of Jh, i.e., Jh > 5 
eV. Since the d — p exchange model can well explain the 
ferromagnetism, it would not be neccessary to have such 
an effective DE model for describing the physics of man- 
ganites. We believe that the d—p exchange -|- (bi)polaron 
model is sufficient to account for the essential physics 
in manganites. 
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